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Abstract. The Early Data Release of the Sloan Digital Sky Survey (SDSS) contains 150 Narrow-Line Seyfert 1 (NLSl) 
galaxies, most of them previously unknown. We present here the study of the X-ray emission from two of these active galaxies 
(SDSS J030639.57 + 000343.2 and SDSS J141519.50 - 003021.6), based upon XMM-Newton observations. The spectral and 
timing characteristics of the two sources are presented and compared against the typical properties of known NLSl galaxies. 
We found that these two NLS 1 are within the dispersion range of the typical values of this class of AGN, although with some 
interesting features that deserve further studies. 
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1. Introduction 

The first studies of Markarian 359 (Davidson & Kinman 1978) 
and of similar objects led to the identification of a pecu- 
liar subclass of Seyfert galaxies: the NaiTow-Line Seyfert 1 
galaxies (NLSl) (Osterbrock & Pogge 1985, Goodrich 1989). 
Such objects presented an intriguing mixture of physical phe- 
nomena, which made them very attractive and popular dur- 
ing the last decade. Their distinctive optical characteristics are 
FWHMinp) < 2000 km/s, weak [OIII] and strong Fe II rel- 
ative to HyS, particularly [OIII]/H^ < 3 (see Grupe 2000, and 
Pogge 2000 for a review). 

Stephens (1989) suggested that X-ray selection could be 
an eflicient method to identify NLSl. Since then, many ef- 
forts have been devoted in analysing X-ray data of NLSl, 
mainly thanks to the satellites ROSAT, ASCA, and more re- 
cently XMM-Newton and Chandra. NLSl typical character- 
istics in the X-ray domain are the presence of a soft X-ray 
excess, probably the high-energy tail of the radiation emitted 
from the accretion disk, a pronounced variability (NLS 1 are the 
most variable AGNs known after blazars), and a steep photon 
index (e.g. Puchnarewicz et al. 1992, Boiler et al. 1996, Wang 
et al. 1996, Grupe et al. 1998, 1999, Leighly 1999a,b). From 
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all these studies, it resulted that NLSl could have the small- 
est central BH mass among all AGN (as low as 10^ Mq) and 
super-Eddington accretion rates (e.g. Wang et al. 1996). 

Since most of NLSl are X-ray selected, it is not clear if 
the above mentioned characteristics are typical of these objects 
or are the result of a biased selection. For example, it is worth 
mentioning that some NLS 1 could have X-ray behaviour sim- 
ilar to Seyfert 1 (see Fig. 8 of Boiler et al. 1996). Therefore, 
it is necessary to try the inverse operation of selection, i.e. to 
study the X-ray characteristics of a large sample of optically se- 
lected NLSl. In this respect, the 150 NLSl found in the Early 
Release of the Sloan Digital Sky Survey (SDSS') by Williams 
et al. (2000) represent a valuable starting point. 

In order to pursue this goal, the XMM-Newton public 
archive has been searched for coincidental position with SDSS 
galaxies of the sample of Williams et al. (2000). Two SDSS 
sources, SDSS J141519.50 - 003021.6 (z = 0.135) and SDSS 
J030639.57 + 000343.2 (z = 0.107), turned out to be in the 
XMM-Newton observed fields. The distance of the sources from 
the boresight was 603" and 170", respectively. For the sake of 
simplicity, they will be thereafter called SDl and SD2, respec- 
tively. These sources were already observed in X-rays only by 
ROSAT (see Williams et al. 2002), and therefore the analysis 
of XMM-Newton data presented here represent the first 2-10 
keV study. 
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Fig. 1. EPIC-PN comparison of source and background spectra for SDl {left) and SD2 {right). The upper data refer to the source 
and background, while the lower one is the background only extracted from a region equal to that of the source. 
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2. XMM - Newton observation and data reduction 

For the processing, screening, and analysis of the data from the 
EPIC MOSl and MOS2 (Turner et al. 2001) and PN cameras 
(Struder et al. 2001), standard tools have been used (XMM SAS 
V . 5 . 4 . 1 and HEAsoft Xspec 11.3.1 and Xronos 5.19) 
and standard procedures described in Snowden et al. (2002) 
followed. The observations of both sources were affected by 
solar soft-proton flares, so that a cleaning was necessary. 

To study variability, the lightcurves from EPIC-PN data 
were extracted, since this detector has the best time resolution 
in full frame mode (~ 73 ms) and the highest throughput. The 
source region was centered in the optical position of the two 
NLSl, a = 14 : 15 : 19.50, 5 = -00 : 30 : 21.6 for SDl, and 
a = 03 : 06 : 39.57, ^ = H-00 : 03 : 43.2 for SD2 (J2000, 
SDSS uncertainty 0.1", Pier et al. 2003), with a radius of 40" 
and 35", respectively. The background to be subtracted in the 
analysis was derived from an annular region with maximum ra- 
dius of 2' for SDl, and from a circular region 1' wide near the 
source for SD2. The different radii of extraction were due to the 
position of the source in the detector chip (closeness to gaps). 

Data from the same regions were also used for the spec- 
tral analysis. The spectra were rebinned so that each energy 
bin contained a minimum of 25 counts, and fitted only in the 
0.5 - 10 keV energy range because of the uncertainties in the 
MOS cameras calibration at lower energies (cf. Kirsch 2003). 
The photon redistribution matrix and the related ancillary file 
were created appropriately with the rmfgen and arfgen tasks 
ofXMM-SAS. 

Despite the high background, and the low statistics (partic- 
ularly for SDl), the signal is significantly higher when com- 
pared to the background level, up to high energy (Fig. [Q. 
Therefore, any feature at high energy can be due to the source. 
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Fig. 2. EPIC-PN global fight curve of SDl, after cleaning for 
soft-proton flares, binned at 100 s. EiTorbars are at Icr. 



3. SDSS J141519.50 - 003021.6 (SD1) 

The source SDl was in the field of view of the ObsID 
0145480101 during the observation performed on 2003 - 08 - 
02. All the EPIC detectors were set in full frame mode and col- 
lected data for an elapsed time of 23567 s. However, the pres- 
ence of soft-proton flares reduced the effective exposure to 6.7 
ks for both MOS and 6. 1 ks for PN. 

The source displays a hint of variability (Fig. with the 
X' probability of constancy less than 7%, and of a 476 s peri- 
odicity (2cr) in the energy band 0.5 - 2 keV. However, the low 
statistics prevents more detailed studies. Longer and unflared 
observations are needed to assess this feature. 

We started the spectral analysis by fitting the average spec- 
trum integrated over the whole efi'ective exposure, with a sim- 
ple power law model and the Galactic absorption along the 
line of sight (A^h = 3.2 x 10^" cm"^. Dickey & Lockman 
1990). This resulted in a steep (F = 2.9 + 0.2) spectrum, but 
with a low statistical quality {x^ = 81.4, dof= 46). The ratio 
data/model clearly shows excesses below 0.6 keV and above 
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Fig. 3. {left) Ratio data/model of the EPIC data of SDl fitted with only a power law absorbed by the Galactic column. Circles are 
MOSl data, crosses are MOS2 data, and squares are PN data. The soft and hard excesses are clearly visible. Data are rebinned 
for clarity, {right) Unfolded spectrum of SDl with the best fit model (partial covering). 



2 keV (Fig. |3}- The excess at high energy cannot be due to 
residual background, as shown in Fig.^ 

The addition of a redshifted thermal component (zbb model 
in xspec) with kT - 0.14+0.02 keV improved significantly the 
fit {x^ - 54.6, dof= 44). In this case, the photon index of the 
power law model decreases to F = 1.8+ 0.2, but still with some 
residuals at high energy. The fit with a gaussian emission line 
fixed at £■ = 6.5 keV with a width cr - I keV was rejected {x^ = 
54.8, dof= 43), and also other models, like the constant density 
ionized model by Ross & Fabian (1993), or the reflection by 
ionized material by Magdziarz & Zdziarski (1995), were not 
successfully constrained. 

A partial covering (pcfabs model in xspec) appears to be 



the best fit model {y^ 
10^2 cm 

r= 3.4: 

0.5 - 10 keV is 4.6 x 10"'^ erg cm^^ s"' 
intrinsic luminosity of 1.2 x lO**"^ erg/s. 



51.0, dof= 44), with A^h = (4!^) x 
^, a covering factor of 0.84 + 0.01, and a photon index 
: 0.2. In this case, the observed flux in the energy band 



corresponding to an 




start Time 12682 3:45:41: 52 Stop Time 12662 14: 2:21: 52 

Fig. 4. EPIC-PN global light curve of SD2, after cleaning for 
soft-proton flares, binned at 1000 s. Error bars are at Icr. 



4. SDSS J030639.57 + 000343.2 (SD2) 

The source SD2 was observed on 2003 - 02 - 1 1 in the field of 
view of the ObsID 0142610101, with an elapsed time of 73918 
s. All the EPIC camera detectors were set in full frame mode. 
Because of the presence of soft-proton flares, the effective ex- 
posures were about 52 ks for both MOS and 38 ks for PN. 

The source does not show statistically significant variabil- 
ity (Fig. 13, the;^^^ probability of constancy being greater than 
92%, depending on the bin width (100 - 1000 s). This appears 
to be unusual, given the well known variability of the NLSl 
in the X-ray energy band: Leighly (1999a) reported variability 
with at least 99% confidence level. The steady flux is, however, 
understandable if the cause of variability is mainly a flare-like 
phenomenon. In this case, we were observing a period without 
flares. 

Therefore, it is reasonable to study only the averaged spec- 
trum integrated over the whole period. The first fit was ob- 
tained with a simple power law model (F = 1 .99 + 0.02), with 



the Galactic absorption toward SD2 {Nh - 6.94 x 10^" cm"^. 
Dickey & Lockman 1990). The result is shown in Fig. |5] soft 
and hard excesses are clearly visible {x^ - 806.9, dof= 561). 
To fit the soft excess we adotped the redshifted black body (zbb 
model in xspec) with kT = 0.083 + 0.004 keV. The photon 
index decreases to F = 1.91 + 0.03, and the fit improves to 
X' - 739.9 for dof= 559. The thermal component is required 
at > 99.99%. 

An excess at energies greater than ~ 5 keV remained in the 
data still after having used diff'erent procedures in the clean- 
ing for soft-proton flares, but the comparison with the back- 
ground (Fig.[0 confirms that the excess is real. To model this 
hard X-ray excess, we tried with a simple large gaussian line 
at £■ = 6.6+°^ keV, cr = 0.8+°^ keV, and equivalent width 
5 10;^234*' This line is required at 99. 15% level, and improves 
the fit to x^ - 724.5 for 556 degrees of freedom. The thermal 
component is now kT — 0.079^QQgj keV and the photon index 
is F = 1 .95;^Q Q^. This is the best fit model for this source (Fig.|5l 
right). The observed flux in the energy band 0.5 - 10 keV is 
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Fig. 5. (left) Ratio data/model of the EPIC data of SD2 fitted with only a power law absorbed by the Galactic column. Circles are 
MOSl data, crosses are MOS2 data, and squares are PN data. The soft and hard excesses are clearly visible. Data are rebinned 
for clarity, (right) Unfolded spectrum of SD2 with the best fit model (power law plus black body and gaussian emission line). 



6.9 X 10 erg cm ^ s \ corresponding to an intrinsic lumi- 
nosity of 2.3 X lO''^ erg/s. 

We tried also a partial covering, obtaining A^h - (19^J^) X 
10^^ cm"^, a fraction of 0.37^^^^, and a photon index T = 
2.07 ± 0.04. The fit, with = 766.7 for 559 degrees of free- 
dom, is worse than the previous one. 

Also the Ross & Fabian (1995) model does not provide a 
good fit (x^ = 749.0, dof= 559). In this case, we found that 
the ionization parameter is log^ = 3.25 ^g gj, with a reflection 
fraction of O.S^'^ l, and T = 1.85 ± 0.04. 

5. Observations at other wavelengths 

5.1. SDSSJ141519.50 -003021.6 (SD1) 

The optical counterparts of SDl have been found in 2MASS 
(Cutri et al. 2003) and USNO Bl (Monet et al. 2003) catalogs. 
The 2MASS position is /?A = 14 : 15 : 19.51 and Dec = -00 : 
30 : 21.1 (J2000, uncertainty 2"), consistent with that of SDSS. 
The source was detected with magnitudes / = 16.7 ± 0.2, H = 
16.5 + 0.2, and K = 15.4 + 0.2. The USNO B 1 catalog provides 
the magnitudes in two periods: (1) Palomar Observatory Sky 
Survey (POSS) I, (1949 - 1965), with emulsion sensitive at 
wavelengths in the range 620 - 670 nm; (2) POSS II (1985 - 
2000), sensitive at 385 - 540 nm. These are: Bi = 19.3 ± 0.3, 
Bl = 19.6 ± 0.3, and Ri = 18.6 + 0.3, R2 = 18.6 ± 0.3, plus 
/ = 17.9 ± 0.3. The detections in the two periods do not show 
variability within the measurement errors. Also in this case, the 
position found (RA = 14 : 15 : 19.5, Dec = -00 : 30 : 21.4, 
J2000) matches that of the SDSS. 

5.2. SDSS J030639.57 + 000343.2 (SD2) 

For this source, in addition to the counterparts in the 2MASS 
and USNO Bl catalog, a radio counterpart was also found in 
the FIRST^ and NVSS^ surveys, both at 1.4 GHz (20 cm). The 

^ http://sundog.stsci.edu/ 

^ http://www.cv.nrao.edu/nvss/ 



FIRST survey was performed between 1993 and 2002, while 
the NVSS was done between 1993 and 1997. The source ap- 
pears to be compact, with a peak flux is 4.4 + 0.4 mJy (NVSS). 

This galaxy appears to be very bright in the optical/infrared 
wavelengths: the 2MASS catalog reports a clear detection with 
allfilters: / = 15.08+0.08, // = 14.12+0.08, = 13.33±0.05. 
The USNO B 1 catalog gives fliese values: Bi = 16.5 ±0.3, B2 = 
14.5 ± 0.3, Ri = 14.7 ± 0.3, /?2 = 13.5 ± 0.3, and 7 = 14.8 ± 0.3, 
showing a substantial degree of variability. Specifically, it is 
worth noting the change of about 2 magnitudes in the B band 
between the two reference periods. 

6. Discussion 

We presented here the first analysis in a wide X-ray en- 
ergy band (0.5 - 10 keV) performed to date on two Narrow- 
Line Seyfert 1 Galaxies from the optically selected sample of 
Williams et al. (2002). By comparing the spectral parameters 
obtained for the present sources with the average values found 
by Leighly (1999a,b), it is clear that the X-ray characteristics 
are in the ranges of the NLS 1 . The mean value of the photon 
index found by Leighly (1999b) is 2.19 ± 0.10, with a disper- 
sion of 0.30!jJ j;^, for NLSl, and 1.78 ± 0.11, with a dispersion 
of 0.29^g for Seyfert 1 active nuclei. These values should be 
compared with the photon index of 3.4 for SDl and 1.95 for 
SD2. SDl appears steeper than the mean value of Leighly, but 
it is worth noting that a steep photon index can stiU be likely 
(see Fig. 8 by BoUer et al. 1996). 

Only SDl displays a certain degree variability, although 
the limited eff^ective exposure (6 ks) prevents further stud- 
ies. SD2 does not show variability, but this could be ex- 
plained with a period of quiescence. To compare the present 
fluxes with the ROSAT observations, we converted, with 
WebPIMMS"^, the PSPC count rates into flie flux in flie band 
0.5 - 2.4 keV, by using the Galactic A'h and the photon in- 
dex r reported by Williams et al. (2002). The calculated ob- 
served flux is (1.4 + 0.4) X lO"'-' erg cm"^ s"' for SDl and 

http : //heasarc . gsf c. nasa . gov/Tools/wSpimms . html 
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(1.1 + 0.2) X 10 '2 erg cm ^ s ' for SD2. The correspond- 
ing fluxes from the present XMM-Newton observations are 
(3+^) X 10"'^ erg cm^^ and (3.2 + 0.1) x 10"'^ erg cm^^ s^^, 
respectively, confirming the variability of both sources, partic- 
ularly for SD2. 

Both sources present an excess at high energy, that could be 
explained by a partial covering model for SDl and by a gaus- 
sian line for SD2. The partial covering has been successfully 
used to explain similar hard X-ray excesses in IH 0707 - 495 
(Grupe et al. 2004, Gallo et al. 2004a, Tanaka et al. 2004). 
However, it is worth stressing that for SDl the statistics is not 
enough to reach firm conclusions. 

The large gaussian emission line used to fit the excess 
in SD2 can be the indication of light bending around a Ken- 
black hole, like, for example, the case of PHL 1092 (Gallo et 
al. 2004b) or reflection from ionized disc, used also for IH 
0707 - 495 (Fabian et al. 2004). Indeed, for SD2, the con- 
stant density ionized disc model by Ross & Fabian (1993) gives 
acceptable fits, even though worse with respect to the simpler 
phenomenological model of the gaussian emission Une. 

Dedicated observations with higher statistics are needed to 
assess the spectral parameters of SDl and SD2. Moreover, in 
order to understand if the characteristisc of NLSl are biased 
by the X-ray selection, it is necessary to enlarge the sample of 
observed galaxies in X-ray from the optically selected sample 
ofWilUamsetal. (2002). 
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